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Extending a Pade approximant method used for studying compactons in the Rosenau-Hyman 
(RH) equation, we study the numerical stability of single compactons of the Cooper-Shepard-Sodano 
(CSS) equation and their pairwise interactions. The CSS equation has a conserved Hamiltonian 
which has allowed various approaches for studying analytically the nonlinear stability of the solu- 
tions. We study three different compacton solutions and find they are numerically stable. Similar 
to the collisions between RH compactons, the CSS compactons reemerge with same coherent shape 
when scattered. The time evolution of the small-amplitude ripple resulting after scattering depends 
on the values of the parameters I and p characterizing the corresponding CSS equation. The 
simulation of the CSS compacton scattering requires a much smaller artificial viscosity to obtain 
numerical stability than in the case of RH compacton propagation. 
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I. INTRODUCTION 

Following their discovery [ ], compactons, or solitary 
waves defined on a compact support, have found diverse 
applications in physics [2, 3], ocean dynamics [4], magma 
dynamics [5, 6], mathematical physics [' -9], nonlinear 
lattice dynamics [10-16], and medicine [17]. Multidimen- 
sional compactons have also been discussed in [18, 19], 
and compact structures have been studied in the context 
of the discrete Burridge-Knopoff model [20], and in the 
context of discrete and continuous Klein- Gordon mod- 
els [21-23]. A recent review of nonlinear evolution equa- 
tions with cosine/sine compacton solutions can be found 
in Ref. 24. 

The K(l,p) compactons discussed first by Rosenau and 
Hyman (RH) are examples of a class of traveling-wave 
solutions with compact support resulting from the bal- 
ance of both nonlinearity and nonlinear dispersion. RH 
discovered these compactons in their studies of pattern 
formation in liquid drops using a family of fully nonlinear 
Korteweg-de Vries (KdV) equations [1], 



u t + {u l ) x + (u p ) 3 



= , 



(1.1) 



where u = u(x, t) is the wave amplitude, x is the spatial 
coordinate and t is time. Equation (1.1) is known as the 
K(l,p) compacton equation. The RH compactons have 
the remarkable soliton property that after colliding with 
other compactons they reemerge with the same coherent 
shape. The collision site is marked by the creation of 
a compact ripple. The positive- and negative-amplitude 
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parts of the ripple decay slowly into low-amplitude com- 
pactons and anti-compactons, respectively [ ]. De Frutos 
et al. showed [25], and Rus and Villatoro confirmed re- 
cently [26], that shocks are generated during compacton 
collisions. 

In general, Eq. (1.1) does not exhibit the usual en- 
ergy conservation law. Therefore, Cooper, Shepard and 
Sodano (CSS) proposed a different generalization of the 
KdV equation based on the first-order Lagrangian [27] 



L(l,p) = 



l^t+jj^-a(4> x ) p (<l> xx ) 2 ]dx, (1.2) 



which leads to the equation: 

u t + u t - 2 u x -p[u p - 1 {u x ) 2 } x + 2a[u p u x } xx =0. (1.3) 

Here, we have u = <j) x . Since then, various other La- 
grangian generalizations of the KdV equation have been 
considered [28-32]. The equation for the solitary waves 
is obtained by substituting u(x, t) = fix — ct) = f(y) 
into Eq. (1.3) and then integrating twice and setting the 
integration constants to zero. One obtains: 



-f 



1(1-1) 



+ a(f) 2 f p = 



(1.4) 



Anti-compacton solutions correspond to the transforma- 
tion / — > — /. Therefore, from Eq. (1.4), we find that for 
anti-compactons to exist, I — p must be an even integer. 
Moreover, when p is odd c changes sign and the anti- 
compacton travels with negative velocity, whereas for p 
even the velocity of the anti-compacton is positive. 

Compacton solutions are constructed by patching a 
compact portion of a periodic solution that is zero at 
both ends to a solution that vanishes outside the com- 
pact region to give a weak solution to the equation. We 
see that for there to be a solution of that type, p < 2 
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and I > p. The condition for a weak solution is that the 
jump across the boundary of the equation of motion at 
Xq where f[xo]=0 is zero. That is 



Di SC [(f) 2 n 



. 



(1.5) 



This is always satisfied if there is no infinite jump in the 
derivative of the function. The stability analysis of the 
solutions relies on the fact that the equation of motion 
for f(y) can be obtained from an Action functional: 



m = J 



dy 



J 2 



f 



a(f') 2 f P 



2 J 1(1-1) 
We recognize this functional as the value of 
P[f]c + H[f], 



(1.6) 



(1.7) 



where P[f] and H[f] are the values of the conserved mo- 
mentum and Hamiltonian respectively for the solitary 
wave /. The once integrated equation of motion for the 
solitary wave is obtained from the equation 



<5$ 

^7 



(1.8) 



and the unintegrated equation of motion for f[y] is given 
by 



S1> 







(1.9) 



For the RH equation, stability of the compacton has 
been demonstrated numerically as well as by a linear sta- 
bility analysis of the radiation induced by the numerical 
method [33]. For the CSS equation, because of the exis- 
tence of a Hamiltonian formulation, various other meth- 
ods of studying nonlinear stability have been explored 
such as Lyapunov stability [30, 34] and stability of the 
solutions under scale transformations [30, 35]. However, 
apart from a numerical study of the evolution and scat- 
tering of the compactons in the generalized CSS equation 
by Cooper, Khare and Hyman [29] using pseudospectral 
methods, there has been no systematic study until now 
of the stability of the compacton solutions to the CSS 
equation. Nor has there been any study of whether the 
solutions that arise from a Hamiltonian dynamical sys- 
tem behave differently from those obeying the four con- 
servation laws of the RH equation [1] (without energy 
conservation). It is this gap in our knowledge that we 
hope to fill by this study. 

To study stability under scale transformations it is suf- 
ficient to study the change in the Hamiltonian for fixed 
momentum P [30]. That is we let 

f{x) -► P l ' 2 f(px), (1.10) 
which leaves P = J dxf 2 /2 unchanged. The Hamiltonian 

1 



H 



dx 



1(1-1 



-f 



H 1 -H 2 , (1.11) 



is then transformed into 

H X (fi) = /35('+3p-2) i J 1 _ pW-2) H2 

The exact solution satisfies: 



dH 



13=1 



= 0. 



This yields 



(I - 2)H 2 = (p + 4)H X 



(1.12) 



(1.13) 



(1.14) 



The second derivative at /3 — 1 can then be written as 



d 2 H 



dp 2 



= ^(p+i^p-i+mi 



Since Hi and H 2 are positive definite we find that the 
solutions are stable to a small scale transformation when 



2 < I <p + 6. 



(1.15) 



This includes all the solutions we will be studying here. 

In a recent paper [36], we performed a systematic 
derivation of a Pade approximants method [37] for cal- 
culating derivatives of smooth functions on a uniform 
grid by deriving higher-order approximations using tra- 
ditional finite-differences formulas. Our derivation con- 
tained as special cases the Pade approximants first intro- 
duced by Rus and Villatoro [26, 33, 38]. We note that 
the L(l.p) compactons feature higher-order nonlinearities 
and terms with mixed-derivatives that are not present in 
the K(p,p) equations. Therefore, in this paper we will 
extend our earlier approach [36], so that we can study 
the compactons that occur in the CSS equation. This 
approach can also be applied to the recent VT general- 
izations of that equation [32] . 

This paper is outlined as follows. In Sec. II, we review 
briefly the main findings with respect to the numerical 
schemes based on Pade approximants derived in Ref. 36. 
Our numerical approach to solving the CSS equation is 
described in Sec. III. In section IV we study numerically 
the stability of several compacton solutions of the CSS 
equation, and we also study the pairwise interactions of 
these compactons. We compare our results on stability 
with our previous numerical study of the if (2, 2) equa- 
tion [36]. We summarize our main findings in Sec. V. 



II. PADE APPROXIMANTS 

We consider a smooth function u(x), defined on the 
interval x G [0, L], and discretized on a uniform grid, 
x m = mh, with m = 0, 1, • • • , M, and h = L/M. Pade 
approximants of order k of the derivatives of u(x) are 
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defined as rational approximations of the form 



,W - 



,(«) - 



,(*») 



(to) _ 



A[E) 
HE) 
B(E) 
HE) 
C(E) 
HE) 
V(E) 
HE) 





\-0(Ax k ) , 


(2.1) 


u m 


\-0(Ax k ) , 


(2.2) 


v> m 


VO{Ax k ) , 


(2.3) 


u m 


\-0{Ax k ) , 


(2.4) 



where we have introduced the shift operator, E, as 

E k u m = u m+k . (2.5) 

Even- and odd-order derivatives require approximants 
that are symmetric and antisymmetric in E, respec- 
tively. The familiar second-order accurate approxima- 
tion of derivatives based on finite-differences are trivial 
examples of Pade approximants 



ME) 
Cx{E) 
Vi{E) 



1 



2Ax 
1 

Ax 2 
1 



2Ax z 
1 

Ax^ 



E- E- 1 
E-2 + E- 1 

E 2 -2E + 2E- 1 - E~ 2 
E 2 -AE + Q-AE- 1 + E~ 2 



(2.6) 
(2.7) 
(2.8) 
(2.9) 



corresponding to Fi{E) = 1 [33]. Still keeping F\{E) = 
1, but incorporating the additional grid points, {x rn ±2\, 
we can obtain fourth-order accurate approximation for 
the derivatives Um and u m , as 



ME) 

He) 



\2Ax 
1 

2Ax 2 



E 2 -8E + 8E- 1 - E~ 



, (2.10) 



E 2 - 6E + 10 - 6E- 1 + E~ 2 



(2-11) 



Previously [3G], we showed on general grounds that 
the Pade-approximant approach allows one to improve 
the numerical representations of only three of the four 
lowest-order derivatives of u(x), when involving only the 
grid points {x m , x m ±i, x m ±2}- To obtain a fourth-order 
accurate approximation of the derivatives, we can ei- 
ther begin by improving the third-order derivative, u m l ^ , 
or the fourth-order derivative, »m'. Because in the 
compacton-dynamics problem [1, 25, 26, 33, 39-41], the 
fourth-order derivative enters only through the artificial 
viscosity term needed to handle shocks, in Ref. 36 we 
chose to improve the approximation corresponding to 

the third-order derivative, lit"', and focused on obtain- 
ing fourth- or higher-order accurate Pade approximants 
of u m , Um \ and u m , on the subset of grid points, 



Using Eqs. (2.3), we introduced an operator, T(E), 
symmetric in E, as 



HE)u^ = 



1 



(E 2 +E- 2 ) +b(E + E- 1 ) +c 



such that 

HE) = d(E) u m + 0(Ax k ) , 
and showed that for 

a = At , b = r - 4 , c = 2(r + 3) , 
we obtain 

Ci(E) ,„.. rt /l 1\ Ax 4 



,/(*") — : - 

m ~ T(E) " 



m Uo t) 



(ix) 



1 2520 t) 



1\ Ax 6 
tI ~2A 



A 

0(Ax 8 ) 



„(«*) 
(2.12) 

(2.13) 

(2.14) 

(2.15) 



Correspondingly, the Pade approximant of the first-order 
derivative, u m \ is obtained as 



ME) 



u Ax 



HE) 

.„(»«) f A 

m V 105 At 



1 \ Ax 6 



4 

0(Ax s 



(2.16) 



with 



A 2 (E) = 



E 2 + 10E - 1QE- 1 - E~ 2 



(2.17) 



24Ax 

and the Pade approximant of the second-order derivative, 

(«) • 
Um , is 



,(«) 



82(E) 
HE) 



Um—U. 



29 



(vi)f_ 



(via) ( 
m \8A0 



V180 
1\ Ax' 



1_\ A.r 1 



(2.18) 
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+ 0(Ax 8 ) , 



with 
B 2 (E) 



1 

6A2; 2 



E z 



2E-6 + 2E- 1 + E- 2 



(2.19) 



and the Pade approximant of the fourth-order derivative, 

(to) . 
Um 1 IS 

^ = Sir Um+u - l) if + ° (Aa;4) ■ (2 - 20) 

In order to numerically study the stability and dynam- 
ical properties of compactons, we will consider a suite of 
different fourth-order accurate approximation schemes, 
to make sure that results are independent of the pecu- 
liarities of a particular approximation scheme. Therefore, 
just like in Ref. 36, we will consider here several sets of 
approximants that mix fourth-order accurate approxima- 



tions for two of the derivatives Um , u. 1 



,,,.(/,„ , and lit"', with a 



4 



sixth-order accurate Pade approximant for the third one, 
together with an "optimal" fourth-order approximation 
scheme that minimized the extent of the radiation train 
in our previous study of K(2, 2) compactons. 

(6,4,4) scheme: This approximation scheme is an 
extension of the scheme introduced by Sanz-Serna et 
al. [25, 39] using a fourth-order Petrov-Galerkin finite- 
element method, and corresponds to choosing r = 30 in 
Eqs. (2.15) and (2.16). We have 



•^[644] (E) 



1 

120 



E 2 + 26E + 66 + 26E~ 



E- 2 . 
(2.21) 



(4,6,4) scheme: A sixth-order accurate approximation 

for Um , can be obtained by requiring r = 180/7. Then, 
we have 



■^[464] {E) — 



1 
720 



IE 1 + 152E + 402 + 152£T 



7E- 
(2.22) 



(4,4,6) scheme: For r = 60, the coefficient of Ax 4 van- 
ishes in Eq. (2.15) and we obtain a sixth-order accurate 
approximation for Um , as 



•^"[446] (E) 



1 

240 



E z + 56£ + 126 + 56£T 



E~ 



(2.23) 



This scheme is an extension of the scheme introduced 
first by Rus and Villatoro [26, 33]. 

(4,4,4) scheme: Finally, for r = 5, the smallest value 
of t leading to integer positive values of a, b, and c, we 
obtain 



•^"[444] (E) 



1 

20 



E z + E + 16 + E~ 



E~ 



(2.24) 



While only leading to a fourth-order accurate approxima- 
tion scheme, the above choice of r was shown to minimize 
the extent of the radiation train in our previous study of 
K(2, 2) compactons. 



III. NUMERICAL APPROACH 

We will apply the numerical schemes based on the Pade 
approximants discussed above to the case of the equation 



u t - c u x +i]u xxxx + 



1 



/ - 



(3.1) 



ap(uP' 1 ul) x +aq(uPu x - 1 ) xx =0, 



where the subscripts t and x indicate partial derivatives 
with respect to t and x, respectively. Here, u{x,t) is 
time evolved in the moving frame of reference with ve- 
locity Co, and in the presence of an artificial dissipation 
(hyperviscosity) term based on fourth spatial derivative, 
i]d 4 u/dx 4 . The hyperviscosity term is needed to damp 



out explicitly the numerical high-frequency dispersive er- 
rors introduced by the lack of smoothness at the edge of 
the discrete representation of the compacton (see e.g. dis- 
cussion in Ref. [29]). The addition of artificial dissipation 
is responsible for the appearance of tails and compacton 
amplitude loss. In our dynamics simulations we choose 
r\ as small as possible to reduce these numerical artifacts 
without significantly changing the solution to the com- 
pacton problem. We note that in the CSS-compacton 
simulations discussed here, we required a hyperviscosity 
value an order of magnitude smaller than the hypervis- 
cosity used in our previous simulations of K(2, 2) com- 
pactons. Unless otherwise specified, we use r\ = 10~ 5 . 

Setting q = 2 in Eq. (3.1), leads to the case of the CSS 
compacton derived from the Lagrangian (1.2), i.e. 

u t - c u x + rju xxxx + — — j- (u l ~ 1 ) x (3.2) 
2a 



ap(v?- 1 u 2 x ) + 



p+ 1 



= 0, 



whereas for a = (q — 1) _1 and q an even integer, we 
obtain the compacton equation for the 'P7'-symmetric 
case discussed in Ref. 32. Hence, even though in the 
following we focus on the discussion of the properties of 
the CSS compactons, the numerical methods developed 
here apply also to the case of T-'T-symmetric compactons. 

To obtain the spatial numerical discretization of 
Eqs. (3.1) and (3.2), suitable for our fourth-order accu- 
rate Pade-approximant approach, we introduce a uniform 
spatial grid in the interval x € [0, L] by defining the grid 
points x m = mAi, with m = 0, 1, • • • , M and the grid 
spacing Ax = L/M. Then, we have 



= T{E) 



dt 



[c A(E) - rjV(E)]; 



A(E) 



1 



(tl m ) 



1-1 



.1-1 

aqB(E) \(u m )P ({u x }. 



ap(u m ) p 1 ({u x } m ) 



(3.3) 



In Eq (3.3), u m (t) is a numerical approximation to 
u(x mi t), and we assume that u m (t) obeys periodic 
boundary conditions, UAf{t) = u (t). Also in Eq. (3.3), 
we introduced the notation u{x m ,t), and {u x } m to de- 
note a numerical approximation to d x u(x m , t). The latter 
is calculated using Eq. (2.10). The optimal discretiza- 
tion for the study of CSS compactons corresponds to 
Eq. (3.3). As such, for q=2, we have 



= F{E) 



du r . 



dt 
A(E) 
2a 



[coA(E) - r)V(E)}, 



1 



k) 



i-i 



p + i 



C(E) \(u m )P +1 



*p{u m ) p 1 ({u x } m y 



(3.4) 



In order to numerically calculate the dynamics, we dis- 
cretized the time-dependent parts of Eqs. (3.3) and (3.4) 
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FIG. 1. (Color online) Study of the CSS compacton stability. CSS results are compared with those obtained in the K(2, 2) 
case. Here, ui and U2 are CSS compactons with velocity-independent width, corresponding to the case I = p + 2, with p — 1 
and p — 2, respectively, whereas U3 is a CSS compacton with velocity-dependent width, corresponding to the case p = 2, I — 3. 
The numerically-induced radiation train results at time t~75 were obtained using the (6,4,4) scheme described in the text. The 
compactons were propagated in their comoving frames (co = c) with At=0.1 and Ax=0.1, 0.05, and 0.025. In all cases, the 
radiation appears to be a numerical artifact that is suppressed by reducing the grid spacing, Ax. This indicates that indeed 
these compactons are stable. 



in Eqs. (3.1) and (3.2) by implementing the midpoint rule in time, similar to previous studies [24, 36]. The resulting 

approximate equation for Eq. (3.3) is 



= F{E) 



,n+l 



At 



c A(E) - r)V(E) 



,71+1 



apA{E) 



,n+l 



0-1) 
aqB{E) 



,71+1 



l-l 



(3.5) 



+ <f(Wl +1 + {' 



9—1-1 



Here we introduced the notations, — u m {t n ) and w^ 1 " 1 = u m (t n + At), to indicate evaluations at two different 
moments of time. 
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FIG. 2. (Color online) Comparison of compacton stability results as a function of numerical scheme. The CSS-compactons 
stability study is shown to be robust with respect to the chosen fourth-order accurate Pade-approximant numerical scheme. 
CSS results are compared with those obtained in the K(2, 2) case. Here we plot the radiation trains obtained at t=75, by 
propagating the compactons in their comoving frames (co = c) with Ai=0.1 and Ax=0.05. 



For CSS compactons the discretization is 



= F{E) 



At 



ap 



A{E) [( 



,71+1 



A(E) - V V(E)] ( 



,n+l 
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71 + 1 
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(3.6) 



IV. RESULTS AND DISCUSSION 



In the following, we discuss the case of the CSS 
compacton equation, given in the laboratory frame by 
Eq. (1.3) or, in a frame moving with velocity cq, by 
Eq. (3.2), where in the latter we set the hyperviscosity 
to zero, 77=0. We study the properties of the three exact 
compacton solutions described in Ref. 27. The first two of 
these compactons correspond to a class of solutions with 



r 



l = p + 2. The width of these compactons is independent 
of the compacton velocity, c, and the compactons have 
the general form 



u(x, t) 



c(p + l)(p + 2) 



cos? 



V4a(p + l)(p + 2)- 
(4.1) 



where we introduced the notation 

£(x, t) = x - xq - (c - co)* , 



(4.2) 
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with xq the position of the compacton maximum at t = 0. 
For p = 1 and a = |, Eq. (4.1) gives the compacton 
solution 

= 3 c cos 2 ^=£(x,i) , |£(x,t)| < V37T, 

whereas for p = 2 and a = 3 we obtain the compacton 
solution 

1 



u 2 (x,t) = V6c cos -£(x,t) , |£(x,t)| < 3tt . (4.4) 
L6 J 

The third compacton to be discussed next corresponds 
to the values, p = 2 and I = 3, and the width of this 
compacton depends on velocity. Choosing a = \, we 
find 



1 



u 3 {x,t) = 3c--e(x,t), \{(x,t)\<3V2c. (4.5) 



G 

Using Eq. (1.4), one can study the possibility that the 
above compactons, u\, ui and U3, have anti-compacton 
counterparts. We infer that the CSS equation corre- 
sponding to the u\ compactons (p = 1, I = 3) allows 
for anti-compacton counterparts traveling with a neg- 
ative velocity, similar to the RH compactons. The ui 
equation (p = 2, I = 4) allows for compact solutions with 
negative amplitude, but these anti-compactons' have a 
positive velocity and travel in the same direction as U2- 
Finally, the U3 CSS equation (p = 2, I = 3) does not 
allow for anti-compacton solutions. 

We will compare results of simulations for the above 
compacton solutions of the Lagrangian (1.2), with results 
of similar simulations for compacton solutions of the RH 
generalization of the KdV equation, Eq. (1.1): For p re- 
stricted to the interval 1 < p < 3, the K(p,p) equation 
allows for compacton solutions of the form [26, 41, 42] 

u c (x,t) = A-r C o S ^[p^x,t)] , |£0M)|<~7r, (4.6) 



where 



2cp p-1 
A = r , P = — — , 7 = 



p+l 



2p 



p-1 



(4.7) 



performed in the case of the K(2, 2) compacton that are 
known to be stable. 

As shown in Fig. 1, the numerical compactons propa- 
gate with the emission of forward and backward propa- 
gating radiation. If the compactons are numerically sta- 
ble, then the amplitude of this radiation train is sup- 
pressed by reducing the grid spacing, Ax, which shows 
that the radiation train is a numerical artifact. In 
Fig. 1, we illustrate results obtained with the (6,4,4) 
Pade-approximant scheme. Here we chose a snapshot at 
t=75 after propagating the compacton in the absence of 
hyperviscosity (?7=0) with a time step, Ai=0.1, and grid 
spacings, Ax=0.1, 0.05, and 0.025. The amplitude of the 
radiation train is at least 4 orders of magnitude smaller 
than the amplitude of the compacton. Using the grid re- 
fining technique, we can show that indeed the radiation 
is a numerically-induced phenomenon. The noise is sup- 
pressed by reducing the grid spacing, Ax, indicating that 
all studied CSS compacton (4.8) solutions are stable. 

These results are robust with respect to the choice of 
the fourth-order accurate Pade-approximant numerical 
scheme. As shown in Fig. 2, the extent and amplitude 
of the radiation train is a characteristic of the chosen 
numerical scheme, and results for the CSS compactons 
are "identical" with results obtained for the K(2,2) com- 
pactons, albeit for a scaling in the amplitude of the ra- 
diation for a given choice of the time step (Ai=0.1) and 
grid spacing (Ax=0.05). This scaling is indicative of the 
higher nonlinearity of the CSS equation relative to the 
K(2, 2) equation, as observed also when one compares 
the results for the K(2,2) and K(3, 3) equations [33]. 

We note that the origin of the radiation observed in 
the propagation of a compacton was shown previously to 
be of numerical origin in the case of the K(p,p) equation 
by Rus and Villatoro [33], who also showed that this 
self-similarity depends strongly on the time-integration 
method [43]. The induced radiation depicted in Figs. 1 
and 2 is similar to that of Ref. [33] and therefore one 
would expect that the self-similarity of the radiation is 
also a feature of the CSS equation. 



B. Pairwise interaction of CSS compactons 



For illustrative purposes, we will consider the case of the 
K(2,2) equation (p = 2), with the exact compacton so- 
lution 



u c (x,t) 



4c 



cos 



1 



£(x,t) , |£(M)I<2tt. (4. 



A. Study of compacton stability 

To numerically demonstrate the stability of the CSS 
compacton solutions, we performed a numerical study of 
the compacton propagation in the compacton comoving 
frame (c = c), using the Pade approximations discussed 
above, and we compare with results of similar simulations 



In the following we will show that the CSS compactons 
also have the soliton property of remaining intact after 
the collisions. The ripple generated following the reemer- 
gence of the CSS compactons decomposes into com- 
pactons with or without anti-compacton counterparts, 
depending on the values of the I and p parameters in 
the corresponding CSS equation. In this context, it is 
important to recall that the u\ compactons are the only 
CSS compactons that have anti-compacton counterparts 
traveling with a negative velocity, similar to the RH com- 
pactons. The U2 CSS equation allows for anti-compacton 
solutions with negative amplitude, but with positive ve- 
locity, traveling in the same direction as u-i- Finally, the 
113 CSS equation does not allow for anti-compacton so- 
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FIG. 3. (Color online) Collision of two CSS compactons, 
ui(x,t), with ci = 1 and C2 = 2. The width of ui com- 
pactons is independent of the compacton velocity and they 
correspond to the choice of parameters, p = 1 and I — 3. The 
simulation is performed in the comoving frame of reference of 
the first compacton, i.e. Co = Ci, using the (6,4,4) scheme and 
a hyperviscosity, r\ — 2 x 10~ 5 . In the left panels, the collision 
is shown to be inelastic, despite the fact that the compactons 
maintain their coherent shapes after the collision: The first 
compacton (ci=l) is "at rest" before the collision occurs. As 
shown in the left panels, after the collision the centroid of this 
compacton changes position. The right panels depict the early 
development of the ripple created in the collision process. 



lutions. No evidence of shock formation accompanying 
the collision was observed. 

All simulations described next involve collisions be- 
tween two CSS compactons with velocities c\ = 1 and 
C2 = 2. The compactons are propagated in the comoving 
frame of reference of the first compacton, i.e. cq = c\, 
using the (6,4,4) Pade approximant scheme. All simu- 
lations were performed in the presence of an artificial 
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FIG. 4. (Color online) At intermediate times, the ripple cre- 
ated as a result of the collision of the two ui(x,t) CSS com- 
pactons depicted in Fig. 3 leads to the emergence of the first 
compacton at t ~ 400. The decay process is very sluggish, 
similar to the RH case. 
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FIG. 5. (Color online) We study the correlations between 
the ripple and the two reemerging compactons: we stop the 
simulation of the collision process at time t — 400, see plot 
labelled [0] in panel (a). We use this snapshot to initialize 
two additional simulations: [1] a simulations in which we drop 
the large compacton, and [2] a simulation in which we drop 
both compactons. In panel (b), we compare results of the 
three simulations at t = 430. We show that the differences 
between [0] and [1] are lower than the order of magnitude 
of the noise induced by the numerical discretization of the 
problem (compare with the noise depicted in Fig. 1. The 
differences between [1] and [2] are of the order of the machine 
precision errors. The above indicate a lack of correlations 
between the ripple and the reemerging compactons. 



hyperviscosity. Unless otherwise stated, the hyperviscos- 
ity value was rj = 10~ 5 , an order of magnitude less than 
the hyperviscosity used in our previous simulations of the 
K(2,2) compacton collisions [36]. 



9 





0.16 ■ 




0.12 - 








0.08 - 




0.04 - 
- 




o.oo - 




-0.04 - 




0.02 




- 

0.00 - 








-0.02 




0.02 




0.00 - 








-0.02 




0.02 




- 

0.00 - 








-0.02 




0.02 




- 

0.00 - 








-0.02 




0.02 




0.00 - 








-0.02 



t=400 



t=400 




t=800 





t=1600 




t=2 1 00 




40 



80 



120 



1 60 



FIG. 6. (Color online) Study of the ripple decomposition 
dynamics at late times. Here we assumed that the correlations 
between the ripple and the two reemerging compactons are 
negligible, as indicated in Fig. 5. For illustrative purposes, we 
shift the ripple, such that the position of the first emerging 
compacton from the ripple, at t = 400, is kept fixed (upper 
panel). In the snapshots depicted in the lower panels, we 
observe the emergence of anti-compacton in the t — 1200 
graph, and the emergence of a second compacton at t = 2100. 



We consider first the collision between two CSS ui (x, t) 
compactons, see Eq. (4.3), with parameters p = 1 and 
I = 3. The width of the u\ compactons is independent of 
the compacton velocity. 

Case 1: p = 1, I = 3. In Fig. 3, we depict a se- 
ries of snapshots of this collision process. Just like in 
the K(2, 2) compacton case, the collision is shown to be 
inelastic, despite the fact that the compactons maintain 
their coherent shapes after the collision. The first com- 
pacton is "at rest" before the collision occurs. As shown 
in the left panels of Fig. 3, after collision this compacton 
emerges with the centroid located at a new spatial posi- 
tion. The early development of the ripple created as a 
result of the pairwise compacton collision is illustrated in 
the right panels of Fig. 3. 

In Fig. 4 we illustrate the emergence of the first com- 
pacton from the ripple. We note the very sluggish decay 
process, just like in the RH-compacton case [1]. 

To demonstrate the lack of correlations between the 
ripple and the two reemerged compactons after collision, 
we use the result of the simulation at t = 400 [denoted 
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FIG. 7. (Color online) Collision of two U2(x,t) compactons, 
with ci = 1 and C2 = 2. The width of U2 compactons is also 
independent of the compacton velocity and they correspond 
to the choice of parameters, p — 2 and 1 — 4. The simula- 
tion is performed in the comoving frame of reference of the 
first compacton, i.e. Co = ci, using the (6,4,4) scheme and 
a hyperviscosity, rj = 10 _J . Results are similar to the ones 
depicted in Fig. 3. 



as [0] in panel (a) of Fig. 5] , to initialize two additional 
simulations: [1] a simulation in which we drop the large 
compacton, and [2] a simulation in which we drop both 
compactons. In panel (b) of Fig. 5, we compare results of 
the three simulations at t = 430. Here, we note that the 
differences between [0] and [1] are lower than the order 
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FIG. 8. (Color online) Dynamics of the ripple resulting from 
the collision of two Ui(x,t) compactons depicted in Fig. 7. 
We compare the shape of the ripple at different times in the 
time propagation, after the ripple "separated" from the com- 
pactons. In order to compare the shapes of the ripple at 
different times, we forced the ripples to cross the x axis at 
the point indicated in the figure. Therefore the x coordinates 
indicated here are only intended to indicate the spatial extent 
of the ripple. 



same direction. As going to later times in this simula- 
tion was deemed too expensive computationally, we chose 
to terminate it before any compacton or anti-compacton 
emerged from the ripple. 

Case 3: p = 2, I = 3. In Fig. 9, we illustrate the 
dynamics of the ripple created as a result of the collision 
of two CSS us(x, t) compactons. The u 3 compactons cor- 
respond to parameters, p — 2 and I = 3, and their widths 
depend on the compacton velocity. We note that the rip- 
ple "decays" in a suite of compactons, without any anti- 
compacton counterparts, as the CSS equation for p = 2 
and I = 3 does not allow for anti-compacton solutions. 
The amplitude of the ripple in this case is much larger 
than in the case of collisions between RH compactons or 
CSS compactons with compacton velocity-independent 
widths, and this may explain why the dynamics of the 
ripple-decomposition process is much faster in the p = 2 
and I = 3 case. 



of magnitude of the noise induced by the numerical dis- 
cretization of the problem (e.g. compare with the noise 
depicted in Fig. 1). The differences between [1] and [2] 
arc of the order of the machine precision errors. 

Assuming that the correlations between the ripple and 
the two reemerging compactons are negligible, we can 
study the dynamics of the ripple at late times. For il- 
lustrative purposes, we shift the ripple, such that the 
position of the first emerging compacton is kept fixed. In 
Fig. 6, we depict snapshots of the ripple decomposition 
for t < 2100. Here, we note the emergence of a first anti- 
compacton in the t = 1200 graph, and the emergence of 
a second compacton at t — 2100. In the p — 1 and / = 3 
case, the emerging compactons and anti-compactons are 
moving in opposite directions relative to the remaining 
ripple, which is considerably reduced in amplitude. 

Case 2: p = 2, I = 4. Similar to the collision process 
depicted in Fig. 3, in Fig. 7 we present a series of time 
snapshots illustrating the collision of two CSS v,2(x,t) 
compactons. The width of the u-i compactons is also 
independent of the compacton velocity and these com- 
pactons correspond to the choice of parameters p = 2 
and I — 4. The results depicted in Fig. 7 are similar to 
those in Fig. 3, albeit for the differences in the shape of 
the emerging ripple. 

The dynamics of the ui ripple is illustrated in Fig. 8. 
In order to compare the shapes of the ripple at differ- 
ent times after the ripple "separated" from the reemerg- 
ing compactons, in Fig. 8 we plot them such that they 
all cross the x axis at the point indicated in the figure. 
The shape of the ripple is shown to be evolving very 
slowly, likely as a result of the fact that in the «2 case 
compacton and anti-compacton solutions travel in the 
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FIG. 9. (Color online) Dynamics of the ripple created as a 
result of the collision of two CSS uz{x,t) compactons. The 
width of it3 compactons depends on the compacton veloc- 
ity and they correspond to the choice of parameters, p = 2 
and 1 = 3. The simulation is performed in the comoving 
frame of reference of the first compacton, i.e. Co = ci, us- 
ing the (6,4,4) scheme and a hyperviscosity, r\ = 10~ 5 . We 
note that the ripple "decays" in a suite of compactons, with- 
out any anti-compacton counterparts. The dynamics of this 
process is much faster than in the case of collision between 
RH compactons or CSS compactons with compacton velocity- 
independent widths. 
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V. CONCLUSIONS 

To summarize, in this paper we presented a systematic 
study of the stability and dynamical properties of CSS 
compactons. Several numerical schemes based on fourth- 
order Pade approximants have been employed and the 
results were found to be independent of the numerical 
scheme. We find that for the propagation of the CSS 
compactons in time using the implicit midpoint rule leads 
to stable results. The simulation of the CSS compacton 
scattering requires a much smaller artificial viscosity to 
obtain numerical stability, than in the case of RH com- 
pactons propagation. 

Based on our study, we verified numerically the con- 
clusion of stability regarding the CSS compactons first 
derived based on criteria such as Lyapunov stability [34] 
and stability of the solutions under scale transforma- 
tions [35]. 

Just like in the case of RH compactons, the CSS com- 
pactons preserve their coherent shapes after the collision. 



The ripple generated following the reemergence of the 
CSS compactons depends on the values of the param- 
eters I and p characterizing the CSS compactons: For 
a given set of I and p values, the ripple decomposition 
gives rise to compactons and anti-compacton counter- 
parts, depending on the presence and character of the 
anti-compacton solutions allowed by the corresponding 
CSS equation. The decomposition of the ripple is much 
faster for a class of CSS compactons for which the width 
of the compacton depends on its velocity. No evidence of 
shock formation accompanying the collision was observed 
after the collisions between CSS compactons. 
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